ANALYTICAL  STUDIES  OF  SELECTED  JOVIAN 
DECAMETRIC  PHENOMENA 


By 

CARL  NIELS  OLSSON 


A  DISSERTATION    PRESENTED  TO  THE    GRADUATE   COUNCIL   OF 

THE    UNIVERSITY   OF    FLORIDA 

IN    PARTIAL   FULFILLMENT  OF  THE   REQUIREMENTS   FOR    THE 

DEGREE    OF   DOCTOR   OF    PHILOSOPHY 


UNIVERSITY  OF  FLORIDA 
1970 


UNIVERSITY  OF  FLORIDA 


3  1262  08552  6332 


To 


Jean,   Kris,   and  Kay 


ACKNOWLEDGEMENTS 

I  wish  to  thank  the  members  of  my  committee,  Drs .  Alex  G.  Smith, 
Guy  C.  Omer ,  F.  Bradshaw  Wood,  and  Daniel  C.  Swanson.   I  am  especially 
indebted  to  Dr.  Alex  G.  Smith,  the  chairman  of  the  committee,  for  his 
help  and  encouragement  throughout  my  association  with  the  radio  astronomy 
group  . 

Many  others  have  been  of  great  help  through  their  suggestions 
and  discussions.   They  are:   Drs.  Thomas  D.  Carr ,  George  R.  Lebo , 
Henry  I.  Register,  Jorma  J.  Riihimaa,  Izzy  Shever ,  David  J.  Kennedy, 
Messrs.  Robert  J.  Leacock,  Jorge  May,  Jorge  Levy,  Ray  Sisler  and  Mrs. 
Dolores  Krausche. 

I  am  thankful  for  the  help  I  received  from  Mrs.  Jan  Schultz  and 
Mrs .  Judy  Lipof sky  with  my  computer  programs . 

As  all  others  before  me,  I  owe  a  great  deal  to  all  of  the  members 
of  the  Florida  and  Chile  groups  who  have  served  in  the  nightly  observations 

I  am  especially  thankful  to  Mr.  W.  W.  Richardson.   He  always  found 
time  to  help  me  with  the  many  aggravating  computer  problems  which  arose, 
and  he  is  also  responsible  for  the  excellent  figures  used  in  this 
dissertation. 

A  special  thanks  is  due  Dr.  Stanley  S.  Ballard,  Chairman  of  the 
Department  of  Physics  and  Astronomy,  for  his  support,  both  moral  and 
financial,  throughout  this  work. 

Finally  I  wish  to  thank  my  wife,  Jean,  for  her  encouragement 

and  for  typing  this  dissertation. 

iii 


TABLE  OF  CONTENTS 

Page 

ACKNOWLEDGEMENTS iii 

LIST  OF  TABLES vi 

LIST  OF  FIGURES vii 

ABSTRACT xi 

Chapter 

I .   INTRODUCTION 1 

II.   THE  DATA  RESERVOIR 9 

Observatories   9 

Instrumentation  and  Raw  Data 9 

Other  Data 12 

III.   LONGITUDE  DISTRIBUTION  OF  THE  MILLISECOND  PULSES  ....  15 

Types  of  Pulses 15 

Io  Influence 22 

IV.   VARIATIONS  OF  ROTATION  PERIOD  BY  SOURCE  AND  BY  FREQUENCY  39 

Introduction  39 

Method  of  Analysis 40 

Frequency  and  Source  Dependence  of  the  Drift  40 

V.   A  STUDY  OF  THE  ROTATION  PERIOD  BASED  ON  STORM 

COMMENCEMENT  TIMES  59 

Introduction  59 

Duncan's  Hypothesis   60 


Chapter  Page 

V.   A  STUDY  OF  THE  ROTATION  PERIOD  BASED  ON  STORM 
COMMENCEMENT  TIMES  (Continued) 

Analysis  of  Florida  and  Chile  Observations  61 

Storm  Durations 77 

VI.   SATELLITE  INFLUENCES  ON  THE  RADIATION   93 

Hypothetical  Satellites   93 

The  Florida  Analysis  94 

Amalthea 99 

VII.   SUMMARY  AND  CONCLUSIONS 100 

APPENDIX 103 

LIST  OF  REFERENCES 117 

BIOGRAPHICAL  SKETCH  120 


LIST  OF  TABLES 

Table  Page 

1.     ANTENNAS  AT  THE  FLORIDA  AND  MAIPU  OBSERVATORIES 11 

2-     DATA  COLLECTED  ON  SLOW -SPEED  RECORDER 13 

3.     LOCATION  OF  THE  SOURCES 41 

4  .     YEARLY  AVERAGE  DRIFT  RATE  OF  SOURCES 46 

5.     LONGITUDE  DIFFERENCES  IN  SOURCE  LOCATION 

AT  DIFFERENT  FREQUENCIES 54 

6  .     DIFFERENCES  IN  LONGITUDES  OF  THE  SOURCES 

AT  FIXED  FREQUENCY 55 

7  .     STORM  STATISTICS 90 


LIST  OF  FIGURES 


Figure  Page 

1.  Typical  reception  probability  histogram  as  a  function 
of  System  III  (1957.0)  CML,  showing  the  approximate 
locations  of  the  four  sources 4 

2.  Location  of  the  University  of  Florida  Radio 

Observatory  and  Chilean  field  stations  10 

3.  Average  pulse  character  as  a  function  of  System  III 

CML  for  18  MHz  during  19  52 16 

4.  Number  of  30-second  periods  containing  character  2 

and  millisecond  pulses  in  1961.2  at  18  MHz 18 

5.  Number  of  30-second  periods  containing  character  2 

and  millisecond  pulses  in  1963.8  at  18  MHz 19 

6.  Number  of  30-second  periods  containing  character  2 

and  millisecond  pulses  in  1964.8  at  18  MHz 20 

7.  Number  of  30-second  periods  containing  character  2 

and  millisecond  pulses  in  1963.8  at  22.2  MHz 21 

8.  Relative  number  of  intervals  of  18  MHz  type  2  pulses 

as  a  function  of  longitude  for  various  positions  of  Io.  .    24 

9.  Relative  number  of  intervals  of  18  MHz  millisecond 
pulses  as  a  function  of  longitude  for  various  positions 

of  Io 26 

10.  Contour  plot  of  18  MHz  character  2  activity  as  a 
function  of  longitude  and  the  position  of  Io  for 

the  1963  Florida -Maipu  apparition 29 

11.  Contour  plot  of  18  MHz  millisecond  activity  as  a 
function  of  longitude  and  the  position  of  Io  for 

the  1963  Florida -Maipu  apparition 31 

12.  Number  of  intervals  including  18  MHz  millisecond 

pulses  versus  the  longitude  of  Io 34 


Figure  Page 

13.  Polar  plot  of  18  MHz  millisecond  pulse  activity 
for  the  indicated  range  of  Io  positions  centered 

on  210° 35 

14.  Polar  plot  of  18  MHz  millisecond  pulse  activity 
for  the  indicated  range  of  Io  positions  centered 

on  270° 36 

15.  Apparent  longitude  drift  of  the  sources  at  15 

and  18  MHz 43 

16.  Apparent  longitude  drift  of  the  sources  at  18 

and  22.2  MHz 44 

17.  Apparent  longitude  drift  of  the  sources  at  18 

and  27.6  MHz 45 

18.  Source  drift  rate  versus  frequency  for  1964.9 48 

19.  Source  drift  rate  versus  frequency  for  1963.9 49 

20.  Drift  rate  of  source  B  as  a  function  of  time 50 

21.  Drift  rate  of  source  A  as  a  function  of  time 51 

22.  Drift  rate  of  source  C  as  a  function  of  time 52 

23.  Positions  of  sources  A,  B,  and  C  at  various 

frequencies  in  1962 57 

24.  Commencement  time  of  storms  histograms  of  source  A 

for  1965.9  to  1968.1  at  18  MHz 63 

25.  Commencement  time  of  storms  histograms  of  source  A 

for  1958.1  to  1964.8  at  18  MHz 64 

26.  Commencement  time  of  storms  histograms  of  source  A 

for  1960.2  to  1967.0  at  18  MHz 65 

27 .  Commencement  time  of  storms  histograms  of  source  A 

for  1958.1  to  1968.1  at  22.2  MHz 66 

28.  Commencement  time  of  storms  histograms  of  source  A 

for  1960.2  to  1967.0  at  22.2  MHz 67 

29.  Merged  histograms  of  commencement  time  of  storms 68 

30.  Histograms  showing  commencement  fronts  with  dashed 

lines  for  1964.8  to  1968.1  at  18  MHz 69 


Figure  Page 

31.  Histograms  showing  commencement  fronts  with  dashed 

lines  for  1958.1  to  1963.8  at  18  MHz 70 

32.  Histograms  showing  commencement  fronts  with  dashed 

lines  for  1964.8  to  1968.1  at  22.2  MHz 71 

33.  Histograms  showing  commencement  fronts  with  dashed 

lines  for  1958.1  to  1963.8  at  22.2  MHz 72 

34.  Location  of  the  commencement  front  of  source  A 

during  ten  apparitions  of  Jupiter 73 

35.  Motion  of  the  center  of  source  A  in  the  storm  commencement 
histograms  for  18  and  22.2  MHz 76 

36.  A  conventional  probability  histogram  78 

37.  Distribution  of  storms  within  source  A  for  19  59.1 

to  1961.6  at  18  MHz 79 

38.  Distribution  of  storms  within  source  A  for  1962.7 

at  18  MHz 80 

39.  Distribution  of  storms  within  source  A  for  1963.8 

at  18  MHz 81 

40.  Distribution  of  storms  within  source  A  for  1964.8 

at  18  MHz 82 

41.  Distribution  of  storms  within  source  A  for  1965.9 

at  18  MHz 83 

42.  Distribution  of  storms  within  source  A  for  1967.0 

at  18  MHz 84 

43.  Distribution  of  storms  within  source  A  for  1959.1 

to  1962.7  at  22.2  MHz 85 

44.  Distribution  of  storms  within  source  A  for  1963.8 

at  22.2  MHz 86 

45.  Distribution  of  storms  within  source  A  for  1964.8 

at  22.2  MHz 87 

46.  Distribution  of  storms  within  source  A  for  1965.9 

to  1967.0  at  22.2  MHz 88 


47.  Over -all  distribution  of  storm  duration  for  source  A  .  . 

48.  A  history  of  mean  storm  duration  


91 


IX 


Figure 
49. 

50, 


Reception  probability  versus  the  position  of 
satellite  A 


Reception  probability  versus  the  position  of 
satellite  A  with  axis  shifted 


51. 

52. 
53. 
54. 
55. 
56. 
57. 
58. 
59. 
60. 
61. 
62. 
63. 
64. 
65. 


Reception  probability  versus  the  position  of 
satellite  A  with  Io-related  data  omitted  .  . 


Florida -Chile  merged  histogram  at  18  MHz  for  1960.2. 
Florida -Chile  merged  histogram  at  18  MHz  for  1961.6. 
Florida -Chile  merged  histogram  at  18  MHz  for  1962.7. 
Florida -Chile  merged  histogram  at  18  MHz  for  1963.8. 
Florida -Chile  merged  histogram  at  18  MHz  for  1964.8. 
Florida -Chile  merged  histogram  at  18  MHz  for  1965.9. 
Florida -Chile  merged  histogram  at  18  MHz  for  1967.0. 
Florida -Chile  merged  histogram  at  22.2  MHz  for  1960.2 
Florida -Chile  merged  histogram  at  22.2  MHz  for  1961.5 
Florida -Chile  merged  histogram  at  22.2  MHz  for  1963.8 
Florida -Chile  merged  histogram  at  22.2  MHz  for  1964.8 
Florida -Chile  merged  histogram  at  22.2  MHz  for  1965.9 
Florida-Chile  merged  histogram  at  22.2  MHz  for  1967.0 
Florida-Chile  merged  histogram  at  27.6  MHz  for  1963.7 


Page 

95 

97 

98 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 


Abstract  of  Dissertation  Presented  to  the  Graduate  Council 

in  Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 

Doctor  of  Philosophy 


ANALYTICAL  STUDIES  OF  SELECTED  JOVIAN 
DECAMETRIC  PHENOMENA 

By 

Carl  Niels  Olsson 

December,  1970 

Chairman:  Alex  G.  Smith 
Major  Department:  Astronomy 

The  Jovian  18  and  22.2  MHz  radiation  data,  which  were  collected 
at  the  University  of  Florida  Radio  Observatory  and  the  Maipu  Radio  - 
astronomical  Observatory  on  high-speed  recorders,  were  separated  into 
L-bursts  and  S-bursts.   These  two  groups  of  data  were  then  independently 
studied  to  determine  their  distributions  in  longitude,  and  also  examine 
to  what  extent  the  bursts  were  influenced  by  the  satellite  Io.   The 
noise  storms  from  sources  B  and  C  contained  a  greater  percentage  of 
S-bursts  than  those  from  source  A.  The  L-bursts  from  source  A  seemed 
to  be  insensitive  to  the  position  of  Io ,  whereas  the  source  A  S-bursts 
showed  an  Io  effect.   Both  types  of  bursts  from  sources  B  and  C  were 
influenced  by  Io . 

All  of  the  Florida  and  Maipu  data  at  the  frequencies  of  15,  18, 
22.2,  and  27.6  MHz  were  investigated  in  an  attempt  to  determine  any 
possible  frequency  or  source  dependence  of  the  sinusodial  drift  of  the 
sources.  The  results  of  this  study  were  negative,  possibly  due  to  the 
difficulty  in  locating  sources  B  and  C  at  the  lower  frequencies  and  for 
apparitions  of  low  activity. 

The  18  and  22.2  MHz  data  were  analyzed  in  order  to  theck  the 


hypotheses  of  Duncan  that  (1)  the  radio  rotation  period  is  different 
from  the  System  III  (1957.0)  value  when  it  is  determined  by  considering 
the  commencement  times  of  source  A  storms  for  each  apparition,  and  (2) 
that  the  drift  is  due  to  a  variation  in  the  duration  of  the  storms. 
From  this  study  it  was  concluded  that  this  method  of  determining  the 
rotation  period  is  questionable  and  that  the  drift  cannot  be  accounted 
for  by  a  variation  in  the  storm  durations. 

Finally,  the  question  of  possible  satellite  influence  on  the 
radiation  by  Amalthea  and  by  two  hypothetical  satellites  was  considered, 
No  significant  influence  was  detected  in  this  study. 


CHAPTER  I 
INTRODUCTION 

The  Jovian  decametric  emission  has  been  the  subject  of  intensive 
study  for  some  fifteen  years,  but  as  yet  no  wholly  acceptable  theory  has 
resulted. 

As  a  preliminary  to  the  study  covered  in  the  following  chapters 
of  this  volume,  one  should  briefly  review  the  present  state  of  our 
knowledge  concerning  the  decametric  radiation  from  Jupiter. 

Although  the  Jovian  radio  signals  are  highly  sporadic,  they  tend 
to  occur  as  noise  storms  of  durations  measured  in  minutes  up  to  very  long 
ones  lasting  several  hours.   Early  attempts  to  correlate  these  storms 
with  System  II  central  meridian  longitude  (hearafter  abbreviated  CML) , 
an  optically  derived  system,  suggested  that  the  storms  occurred  only 
when  certain  longitudes  were  on  the  central  meridian;  that  is,  there 
seemed  to  be  three  or  four  source  regions  associated  with  the  planet. 
Further  study  indicated  that  the  radio  sources  drifted  in  longitude 
in  System  II . 

As  more  data  were  collected  and  studied,  it  became  apparent  to 
many  groups  of  investigators  (Shain  1956;  Gallet  1961;  Burke  19  57; 
Garnder  and  Shain  1958;  Carr  et_  a_l .  1958,  1961;  Smith  and  Carr  1959; 
and  Douglas  1960,  1960b,  1962)  that  the  sources  could  be  used  to 
define  a  radio  rotation  period  --  that  is,  a  longitude  system  that 
would  minimize  any  apparent  drift  of  the  sources.   Based  on  the  above 
studies,  a  System  III  (1957.0)  longitude  with  a  period  of  9  55m29!37 


(slightly  shorter  than  System  II)  was  tentatively  adopted  by  the  Inter- 
national Astronomical  Union  (IAU  1962)  in  1962  for  specifying  the 
CML  of  radio  features  of  Jupiter. 

At  present  there  are  three  strong  sources  and  one  rather  weak 
source  identified;  they  are  shown  in  Figure  1.   Their  names  and  approximate 
locations  are  as  follows: 

Source  A  or  Main  source  between  longitudes  200° -300° 

Source  B  or  Early  source  between  100° -200° 

Source  C  or  Late  source  between  300° -20° 

Source  D  or  Fourth  source  between  20° -100° 
In  1962  (Douglas  &  Smith  1963a  Smith  et_  al_.  1965a)  evidence 
appeared  that  the  sources  were  drifting  in  System  III,  implying  a  rather 
sudden  change  in  the  radio  rotation  period.   The  sources  drifted  steadily 
to  higher  longitudes  at  about  10°/year  until  about  1965,  when  the 
sense  of  drift  reversed.   This  apparent  source  drift  will  be  considered 
in  detail  in  following  chapters  . 

The  decametric  noise  storms  consist  of  relatively  short  bursts 
of  radiation.   Most  of  the  bursts  fall  into  two  categories  —  the 
L -bursts,  by  far  the  most  common,  with  durations  from  a  few  tenths 
of  a  second  to  several  seconds  and  the  more  rapid  S -bursts,  whose 
durations  are  measured  in  milliseconds.   On  rare  occasions  a  third 
type  has  been  observed,  lasting  tens  of  seconds.   It  now  appears  that 
the  fine  structure  of  the  L-bursts  is  due  to  scintillations  produced 
in  the  interplanetary  medium  (Douglas  and  Smith  1967)  ,  whereas  there 
is  strong  evidence  that  the  S -burst  structure  is  generated  at  Jupiter 
(Smith  et_  al_.  1966)  . 

Several  groups  have  studied  the  morphology  of  the  fine  structure 
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of  the  bursts  using  both  fixed -frequency  receivers  and  wide -band  receivers 
These  results  are  reviewed  in  the  recent  paper  of  Carr  and  Gulkis  (1969) 
and  also  in  one  by  Riihimaa  et_  al_.  (1970)  . 

The  longitudinal  distribution  of  the  L-  and  S -bursts  has  been 
studied  by  the  Florida  State  University  group  (Baart,  Barrow,  and  Lee 
1966;  Barrow  and  Baart  1967)  and  the  University  of  Florida  group  (Smith 
et  al.  1966;  Olsson  and  Smith  1966  and  1967)  and  will  be  discussed  in 
detail  in  Chapter  3  of  this  dissertation. 

One  of  the  most  exciting  discoveries  concerning  Jupiter  came  in 
1964,  when  the  Australian  E.  K.  Bigg  (1964)  showed  that  the  inner  most 
Galilean  satellite,  Io,  has  a  very  marked  influence  on  the  probability 
of  receiving  decametric  radiation  from  the  planet.   Bigg  found  that 
the  probability  of  receiving  emissions  from  source  B  was  greatly  enhanced 
when  Io  was  in  a  position  in  its  orbit  of  about  90°  from  geocentric 
superior  conjunction,  and  from  sources  A  and  C  when  Io  was  at  240°. 
The  Io  effect  is  more  pronounced  for  sources  B  and  C  than  for  source  A. 
It  has  also  been  shown  that  the  Io  influence  is  strongest  at  the  higher 
frequencies  (27-30  MHz)  and  becomes  progressively  weaker  at  lower 
frequencies  (Lebo  et_  al.  1965;  Register  1968) . 

Bigg's  discovery  of  the  Io  effect  inspired  others  (Lebo  et  al. 
1965;  Dulk  1967;  Tiainen  1967;  Register  1968)  to  investigate  possible 
influences  on  the  emissions  by  the  three  other  Galilean  satellites. 
However,  to  date  any  effect  they  may  have  is  questionable.  Moon  V  or 
Amalthea,  the  satellite  closest  to  Jupiter,  and  two  as  yet  hypothetical 
moons  (Bigg  1969;  Douglas  and  Bozyan  1970)  have  also  been  studied  for 
possible  influences,  with  negative  results. 

Soon  after  the  discovery  of  Jupiter's  radio  emission  the  co- 


discoverers,  Burke  and  Franklin  (1956),  also  showed  that  the  bursts 
are  most  often  circularly  or  elliptically  polarized  in  the  right-hand 
sense,  which  strongly  suggested  that  Jupiter  is  surrounded  by  a  magne- 
toionic  medium.   The  existence  of  a  Jovian  magnetosphere  was  soon 
confirmed  by  the  many  observers  studying  the  decimeter  radiation  from 
the  planet.   The  analysis  of  the  decimeter  radiation  is  basically 
responsible  for  our  present-day  model  of  the  Jovian  magnetic  field, 
giving  such  information  as  the  orientation  of  the  field  with  respect 
to  the  rotational  axis,  direction  of  the  field,  and  an  approximate 
idea  of  its  strength. 

Polarization  studies  in  the  decameter  range  have  been  carried 
out  by  several  observers  (Carr  et_  al_.  1965b;  Barrow  and  Morrow  1968; 
Sherrill  1965;  Kennedy  1969) .   One  finding  of  particular  importance 
is  that  while  at  the  higher  frequencies  the  polarization  is  predominately 
right-hand,  it  becomes  more  and  more  left-hand  with  a  decrease  in 
frequency.   This  fact  has  been  useful  in  studying  Io -related  and  non- 
Io -related  emission, in  deducing  which  hemisphere  the  emission  comes 
from,  and  even  in  helping  to  locate  the  fourth  or  D  source.  A  complete 
and  up-to-date  review  of  the  polarization  phemonema  can  be  found  in 
the  dissertation  of  Kennedy  (1969)  and  the  paper  by  Carr  and  Gulkis  (1969) 

Recent  refinements  of  technique  in  very  long  baseline  (VLB) 
inter ferome try  (Carr  et_  al_.  1965a;  Slee  and  Higgins  1966;  Brown  et  al. 
1968,  Block  et_  al_.  1970;  Dulk  1970)  have  made  it  possible  to  set  an 
upper  limit  on  the  sizes  of  the  sources.   Dulk  (1970)  has  recently 
quoted  an  upper  limit  for  an  incoherent  source  as  O'.'lO  or  400  km  at 
the  planet,  using  an  interferometer  with  a  baseline  of  487000X  operating 
at  34  MHz.  Hopefully,  VLB  interf erometry  can  also  determine  the  positions 


of  the  sources  relative  to  the  disc  of  the  planet  if  one  can  indepen- 
dently determine  the  location  of  the  central  fringe  of  the  interfero- 
meter pattern. 

Any  acceptable  model  of  Jupiter's  decametric  radiation  must  be 
able  to  account  for  all  of  the  above  observational  results.  The  two 
discoveries  which  have  caused  most  theories  the  greatest  difficulty 
are  the  apparent  drifting  of  the  sources  in  longitude  and  the  Io  effect. 
The  recent  model  of  Goldreich  and  Lynden-Bell  (1969)  fits  most  of  the 
data  and  seems  to  be  the  most  promising  at  present. 

Of  the  above-mentioned  observations  and  analyses,  the  present 
study  is  concentrated  in  the  three  broad  areas  of  the  longitudinal 
distribution  of  the  pulses,  source  drifts,  and  satellite  influences. 
In  particular,  the  following  will  be  considered: 

1.  The  longitude  distribution  of  L-  and  S -bursts 

2.  Influence  of  Io  and  L-  and  S -bursts 

3 .  Frequency  dependence  of  source  drift 

4.  Variation  of  drift  among  sources 

5.  Duncan's  (1967)  hypothesis  of  source  drift 

6.  Durations  of  the  storms 

7.  Influence  of  Amalthea  and  the  two  "new"  satellites 
hypothesized  by  Bigg. 

The  investigations  listed  above,  and  to  be  described  in  later 

chapters,  have  led  to  several  publications  in  the  open  literature.   In 

particular,  item  1  is  summarized  in  the  following  papers: 

"Jovian  Rotation  Periods  and  the  Origin  of  the  Decametric 
Burst  Structure,"  A.  G  .  Smith,  G.  R.  Lebo ,  C.N.  Olsson, 
W.  F.  Block,  N.  F.  Six,  and  T.  D.  Carr ,  Proceedings  of 
CALTECH-JPL  Lunar  and  Planetary  Conference,  California 
Institute  of  Technology  (1966),  pp.  128-133 


"Decametric  Radio  Pulses  from  Jupiter:   Characteristics," 
C.  N.  Olsson  and  A.  G.  Smith,  Science,  153,  289-290  (1966) 

"Influence  of  the  Terrestrial  Environment  on  the  Temporal 

and  Statistical  Characteristics  of  Jovian  Decametric  Radiation," 

A.  G.  Smith,  W.  F.  Block,  W.  A.  Morton,  G.  R.  Lebo ,  T.  D.  Carr , 

and  C.  N.  Olsson,  Radio  Science,  1  (new  series),  1167-1171 

(1966) 

"Delineation  of  Jovian  Decametric  Sources  by  Millisecond 
Pulses,"  C.  N.  Olsson  and  A.  G.  Smith,  Nature,  214,  999-1001 
(1967). 

Item  2  is  summarized  in  the  paper  "Influence  of  Io  on  Millisecond 

Pulses  of  Jovian  Decametric  Radiation,"  C.  N.  Olsson  and  A.  G.  Smith, 

Bull .  Am.  Phys .  Soc . ,  11,  853  (1966).   Items  3  and  4  are  summarized  in 

the  papers  "Jovian  Rotation  Periods,"  C.N.  Olsson,  A.  G.  Smith,  and 

G.  R.  Lebo,  Bull.  Am.  Phys.  Soc.  ,  11_,  513  (1966)  and  "Apparent  Drift  of 

the  Radio  Rotational  Period  of  Jupiter,"  C.  N.  Olsson  and  A.  G.  Smith, 

Astron.  J.,  73,  S30  (1968).   Items  5  and  6  are  summarized  in  the  following 

papers  : 

"Analysis  of  the  Jovian  Decametric  Radiation  Based  on  the 
Beginning  of  Storms,"  C.  N.  Olsson,  H.  I.  Register,  and 
A.  G.  Smith,  Bull.  Am.  Phys.  Soc,  13_,  1714  (1968) 

"Analysis  of  the  Jovian  Decametric  Source  Drift  Based  on 
the  Beginning  of  Storms,"  C.N.  Olsson,  H.  I.  Register,  and 
A.  G.  S  mith,  Bull.  Am.  Ast .  Soc.  ,  1_,  202  (1969) 

"Commencement  Times  and  Durations  of  Jupiter's  Radio  Noise 
Storms,"  C.N.  Olsson,  A.  G.  Smith,  H.  I.  Register,  and 
J.  May,  Icarus,  U,  212  (1969) 

"Jupiter's  Decametric  Rotation  Period  as  Defined  by  the 
Commencements  of  Source  A  Storms,"  C.  N.  Olsson,  A.  G.  Smith, 
H.  I.  Register,  J.  May,  Bull.  Am.  Phys.  Soc.  ,  15_,  203  (1970). 


CHAPTER  II 

THE  DATA  RESERVOIR 

Observatories 

The  University  of  Florida  maintains  two  major  radio  observatories; 
one  is  located  in  Florida  and  another  in  Chile.   The  Florida  observatory 
was  situated  in  Gainesville  from  December  of  1956  until  the  summer  of 
1967,  when  it  was  moved  to  a  more  remote  site  in  Dixie  County.   In  19  59 
the  University  of  Florida,  with  the  cooperation  of  the  University  of 
Chile,  established  a  southern  hemisphere  radio  observatory  near  Maipu, 
Chile,  a  suburb  of  Santiago.   These  two  observatories,  located  approxi- 
mately the  same  distance  from  the  equator  but  in  opposite  hemispheres, 
assure  the  best  possible  monitoring  of  Jupiter  year  after  year  as  the 
planet  changes  declination.   In  1964  a  field  station  was  established 
in  an  Andean  valley  about  300  miles  north  of  Santiago  near  the  town  of 
Huanta .   The  locations  of  these  observatories  are  indicated  in  Figure  2. 
All  of  the  data  used  in  the  following  analyses  were  collected  at 
either  the  University  of  Florida  Radio  Observatory  or  the  Maipu  Radio  - 
astronomical  Observatory. 

Instrumentation  and  Raw  Data 

The  radio -astronomical  instrumentation  at  the  observatories  includes 
the  antenna  systems  listed  in  Table  1.   The  signal  from  each  antenna  is 
amplified  by  a  high  quality  communications  receiver  and  recorded  on  pen 


10 


Figure  2.  Location  of  the  University  of  Florida  Radio 
Observatory  and  Chilean  field  stations. 
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TABLE  1 
ANTENNAS  AT  THE  FLORIDA  AND  MAIPU  OBSERVATORIES 


Location 

Frequency 

Type 

No.  of  Elements 

Florida 

12.5  MHz 

H  Broadside 

4 

Florida 

15  MHz 

Yagi 

5 

Florida 

18  MHz 

2  Yagis 

5 

Florida 

22.2  MHz 

Yagi  Polarimeter 

5 

Florida 

27.6  MHz 

Yagi 

7 

Maipu 

12  MHz 

Yagi  Polarimeter 

4 

Maipu 

15.8  MHz 

Yagi  Polarimeter 

3 

Maipu 

18  MHz 

Yagi 

5 

Maipu 

22.2  MHz 

Yagi  Polarimeter 

5 

Maipu 

27.6  MHz 

Yagi  Polarimeter 

7 

Maipu 

30  MHz 

Broadside 

15 

Maipu 

32  MHz 

Broadside 

15 

Maipu 

34  MHz 

Broadside 

15 

Maipu 

36  MHz 

Broadside 

15 

Maipu 

38  MHz 

Broadside 

15 

Maipu 

40  MHz 

Broadside 

15 

12 


recorders.   All  of  the  data  except  those  used  in  Chapter  III  were 
recorded  on  Texas  Instruments  Dual  Rectiriter  recording  oscillographs, 
which  are  slow-speed  recorders.   The  data  in  Chapter  III  were  recorded 
on  high-speed  pen  recorders,  Brush  Mark  II  oscillographs.   The  slow- 
speed  records  were  reduced  and  punched  on  computer  cards  using  the 
method  described  by  Register  (1968) .   The  method  used  in  reducing  the 
high-speed  records  is  outlined  in  Chapter  III.   Table  2  lists  information 
concerning  these  data. 

Other  Data 


Throughout  this  dissertation  the  Jovian  longitude  system  used 
is  the  System  III  (1957.0)  as  defined  by  the  International  Astronomical 
Union  (1962).   This  system  rotates  with  a  constant  period  of  9  55  29.37. 
The  periods  of  the  Jovian  satellites  used  in  this  study  were  obtained  from 
the  American  Ephemeris  and  Nautical  Almanac. 
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TABLE  2 
DATA  COLLECTED  ON  SLOW -SPEED  RECORDER 


Frequency    Mean  Epoch 
12  MHz       1966.0 
15  MHz 


18  MHz 


22.2  MHz 


Station    Listening  Hours    Activity  Hours 
Florida  590  102 


176 
49 

102 

107 
45 

109 

42 

13 

15 

28 

72 

162 

150 

206 

182 

160 

284 

191 

50 

176 

46 

124 

13 

11 

6 

11 

31 

33 

59 

75 

110 

143 

103 

36 

81 

28 

57 


1961.5 

Maipu 

934 

1962.7 

Florida 

610 

1963.8 

Florida 

830 

1964.8 

Florida 

1530 

1965.8 

Florida 

540 

1967.0 

Florida 

830 

1957.1 

Florida 

195 

1958.1 

Florida 

278 

1959.1 

Florida 

376 

1960.2 

Florida 

324 

1960.4 

Maipu 

779 

1961.6 

Florida 

1145 

1961.6 

Maipu 

1190 

1962.6 

Maipu 

1200 

1962.7 

Florida 

1115 

1963.7 

Maipu 

1130 

1963.8 

Florida 

1570 

1964.8 

Florida 

1630 

1964.9 

Maipu 

770 

1965.9 

Florida 

1550 

1966.0 

Maipu 

735 

1967.0 

Florida 

1100 

1967.0 

Maipu 

635 

1958.1 

Florida 

366 

1959.2 

Florida 

3  20 

1960.2 

Florida 

290 

1960.4 

Maipu 

964 

1961.5 

Florida 

610 

1961.5 

Maipu 

1190 

1962.8 

Florida 

1030 

1963.7 

Maipu 

1000 

1963.8 

Florida 

1590 

1964.8 

Florida 

1610 

1964.9 

Maipu 

675 

1965.9 

Florida 

1580 

1965.9 

Maipu 

620 

1967.0 

Florida 

1250 

14 


18  MHz 


22.2  MHz 


TABLE  2  (CONTINUED) 


Frequency 

Mean  Epoch 

Station 
Florida 

Listening  Hours 
278 

Activity  Hours 

27.6  MHz 

1959.2 

1 

1960.2 

Florida 

417 

3 

1961.6 

Florida 

960 

25 

1962.5 

Maipu 

500 

76 

1962.7 

Florida 

700 

26 

1963.7 

Maipu 

1140 

32 

1963.8 

Florida 

1520 

39 

1964.8 

Florida 

1410 

24 

1965.9 

Florida 

1370 

21 

1965.9 

Maipu 

880 

15 

1967.0 

Florida 

1180 

14 

DATA  COLLECTED  ON  HIGH-SPEED  RECORDERS 


1961 

.5 

Florida 

1963 

.7 

Maipu 

1963 

.8 

Florida 

1964 

.8 

Florida 

1964 

.9 

Maipu 

1953 

.8 

Florida 

1964 

.9 

Maipu 

CHAPTER  III 

LONGITUDE  DISTRIBUTION  OF  THE  MILLISECOND  PULSES 

Types  of  Pulses 

The  sporadic  decametric  radio  storms  from  Jupiter  are  received 
as  bursts  of  pulses .   Three  distinct  types  of  pulses  are  discernible 
when  one  studies  the  records  from  a  high-speed  recorder  providing 
moderately  good  time  resolution  (speed  of  5mm/sec) .   On  very  rare 
occasions  there  are  pulses  lasting  tens  of  seconds.   Most  of  the 
emission  consists  of  pulses  which  last  from  tenths  of  a  second  to  several 
seconds;  these  are  generally  referred  to  as  "normal"  or  "swishing" 
pulses.   The  third  type  have  durations  measured  in  milliseconds;  because 
of  the  distinctive  sound  produced  by  a  train  of  these  pulses,  they 
have  been  described  as  "popping"  or  "spitting".   For  this  study  a  pulse 
character  number  was  assigned  to  the  three  types  of  pulses.   Character  1 
is  the  very  long  pulses,  2  the  normal  pulses  and  3  the  millisecond 
pulses . 

With  the  help  of  the  high-speed  pen  recordings,  a  character  number 
was  assigned  to  each  30-second  period  of  Jovian  emission.   Periods 
containing  pulses  of  mixed  types  were  given  half -integral  character 
numbers;  if,  for  example,  a  period  contained  pulses  of  types  2  and  3, 
it  was  assigned  a  character  number  of  2.5. 

Figure  3  shows  the  average  character  number  as  a  function  of 
System  III  Jovian  longitude.  For  convenience  in  analysis  the  planet 
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was  divided,  as  usual,  into  5°  zones  of  longitude.   In  determining  the 
average  character  for  a  given  5°  longitude  zone,  a  computer  program 
summed  the  number  of  30 -second  intervals  of  emission  that  fell  within 
that  zone,  weighting  each  by  its  own  character  number,  and  this  sum 
was  then  divided  by  the  total  number  of  intervals.   It  should  be  remem- 
bered that  a  character  above  2.0  means  that  at  least  some  millisecond 
radiation  occurred  for  that  zone.  Thus  the  figure  suggests  that  the 
millisecond  bursts  tended  to  cluster  in  the  longitude  region  from  110° 
to  200°,  or  in  the  general  vicinity  of  decametric  source  B,  the  same 
conclusion  reached  in  the  preliminary  study  of  the  1962  data  (Smith 
et  al.  1966) .   There  is  an  indication  of  further  clustering  near  source  C 
(310°  to  360°) . 

It  is  perhaps  even  more  instructive  to  analyze  the  pulses  in 
the  manner  shown  in  Figures  4  through  7.  Here  the  computer  has  been 
used  to  total  the  number  of  30 -second  intervals  of  each  character 
(1,  1.5,  2,  2.5,  3.0)  that  occurred  while  each  5°  longitude  zone  was 
on  the  central  meridian.   In  the  figures,  types  2.5  and  3  have  been 
combined  to  produce  the  shaded  regions  of  millisecond  radiation,  while 
the  unshaded  histograms  represent  emission  of  type  2  only.  No  corre- 
lation was  attempted  for  types  1  and  1.5,  since  these  classes  formed 
less  than  one  per  cent  of  the  18  MHz  data  analyzed. 

The  striking  feature  of  these  diagrams  is  that  while  the  milli- 
second radiation  is  concentrated  in  the  familiar  sources  A,  B,  and  C, 
it  seems  to  occur  with  roughly  equal  probability  in  all  three  sources. 
On  the  other  hand,  type  2  emission  is  heavily  concentrated  in  source  A; 
in  1963,  for  example,  74%  of  the  character  2  activity  arose  from  this 
source.  This  result  of  course  explains  the  form  of  Figure  3.   Millisecond 
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pulses  represent  about  one -third  of  the  total  activity  for  sources  B 
and  C,  but  less  than  one -sixth  of  the  activity  for  source  A,  and  thus 
the  average  character  is  depressed  in  the  A  region.   It  is  encouraging 
that  Baart,  Lee,  and  Barrow  (1966)  have  reported  results  highly  similar 
to  Figures  4  through  7,  although  their  study  was  based  on  a  much  less 
extensive  body  of  data . 

A  cross -comparison  of  Figures  4  through  7  suggests  that  there 
may  be  a  temporal  variation  in  the  relative  amount  of  millisecond  pulse 
activity.   For  example,  1963,  18%  of  all  of  the  18 -MHz  activity 
intervals  contained  type  3  pulses,  but  in  1964  the  figure  dropped  to 
8.5%.   The  planet^wide  average  probability  of  emission  for  18-MHz 
pulses  of  all  types  declined  from  0.173  in  1963  to  0.125  in  1964. 

It  is  also  interesting  to  note  the  differences  in  Figures  5 
and  7,  which  are  plots  of  the  same  apparition,  1963,  but  at  different 
frequencies,  18  and  22  MHz  respectively.   The  average  probability  of 
receiving  Jovian  emission  at  22  MHz  is  about  467o  that  of  receiving  it 
at  18  MHz  for  this  apparition  (Register  1968) .   It  is  noted  in  comparing 
the  histograms  that  at  18  MHz  a  far  greater  percentage  of  the  total 
emission  contains  millisecond  pulses  than  at  22  MHz;  in  fact,  five  and 
one -half  times  more  millisecond  pulses  were  received  on  18  MHz  than  at 
22  MHz  for  this  apparition. 

Io  Influence 


Since  it  is  now  well  established  that  the  position  of  the  inner- 
most Galilean  satellite,  Io ,  strongly  influences  the  probability  of 
Jovian  decametric  emission  (Bigg  1964) ,  an  effort  was  made  to  see  if 
this  effect  extends  to  pulse  character.   Figures  8  and  9  show  profiles 


Figure  8.  Relative  number  of  intervals  of  18  MHz  type  2 
pulses  as  a  function  of  longitude  for  various 
positions  of  Io . 
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Figure  9.  Relative  number  of  intervals  of  18  MHz  millisecond 
pulses  as  a  function  of  longitude  for  various 
positions  of  Io. 
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of  18 -MHz  activity  versus  the  System  III  longitude  of  the  central 
meridian  for  the  normal  type  2  and  for  millisecond  pulses,  respectively; 
in  each  case  the  profile  is  drawn  for  eight  different  positions  of  Io 
relative  to  the  earth -Jupiter  line.   The  points  for  these  smoothed 
histograms  were  obtained  from  a  merge  of  the  Florida  and  Chile  data 
from  1961  to  1964.   Figure  8  suggests  that  the  type-2  activity  of 
source  A  was  relatively  insensitive  to  the  position  of  Io;  while  it 
shows  the  well-known  enhancement  when  Io  was  roughly  90°  or  240°  from 
superior  geocentric  conjunction,  a  significant  amount  of  activity  was 
always  present  that  was  apparently  unrelated  to  Io.   On  the  other  hand, 
Figure  9  implies  that  the  millisecond  activity  from  source  A  was 
somewhat  more  sensitive  to  the  location  of  the  satellite,  peaking  when 
Io  was  near  western  elongation;  most  striking  of  all  is  the  sharp 
peaking  of  millisecond  activity  in  source  B  when  Io  was  near  90°,  and 
in  source  C  when  the  satellite  was  near  240°.   Another  method  of  dis- 
playing the  marked  effect  Io  has  on  the  millisecond  pulses  versus  the 
normal  ones  is  shown  in  Figures  10  and  11.   These  are  contour  plots 
similar  to  those  used  by  Bigg  (1964)  and  Lebo  et_  al_.  (1965)  in  their 
investigations  of  the  Io  effect.   To  generate  these  plots,  a  matrix 
was  set  up  with  the  seventy -two  5°  intervals  of  Io  from  superior 
conjunction  as  the  ordinate,  and  the  seventy -two  5°  intervals  of  the 
longitude  of  the  central  meridian  as  the  abscissa,  and  the  computer 
then  summed  the  number  of  intervals  of  normal  and  millisecond  activity 
for  each  pair  of  intervals.   The  relative  numbers  are  indicated  by  the 
different  shadings  in  Figures  10  and  11.   By  comparing  the  two  figures 
it  is  once  again  evident  that  the  Io  influence  is  far  greater  for  the 
millisecond  pulses  than  for  the  normal  ones. 
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This  sensitivity  of  sources  B  and  C  to  the  position  of  Io  led 
to  the  study  shown  in  Figure  12,  an  unexpected  result  of  which  was  a 
rather  clear-cut  confirmation  of  the  division  of  source  B  into  two 
subsources,  B,  and  B  ,  an  effect  only  suggested  by  earlier  work  (Smith 
et  al .  1962) .   This  analysis  utilized  only  those  intervals  containing 
millisecond  pulses  at  18  MHz  for  the  1963  apparition.   The  upper  curve 
is  based  on  radiation  received  while  source  B1  was  on  the  central 
meridian,  the  second  curve  on  intervals  when  B  was  on  the  central 
meridian,  and  so  on.   The  abscissa  is  the  System  III  longitude  of 
the  sub-Io  point  on  Jupiter.   It  is  evident  that  B1 ,  B  ,  and  C,  which 
were  sensitive  to  Io's  position  with  respect  to  the  earth -Jupiter  line, 
are  equally  sensitive  to  the  satellite's  System  III  longitude.   On 
the  other  hand,  A,  which  was  insensitive  to  Io's  position  with  respect 
to  the  earth -Jupiter  line,  is  also  relatively  insensitive  to  the 
satellite's  movement  in  System  III.   Notice  that  distinctly  different 
System  III  longitudes  of  Io  stimulate  B..  as  compared  with  B  .   It  is 
interesting  that  Figure  12,  based  entirely  on  millisecond  pulses, 
agrees  with  Dulk's  (1965)  observation  that  normal  Jovian  emission  is 
stimulated  only  when  Io  is  over  the  hemisphere  110°  <  \  <   290°. 

The  angular  relationships  are  more  striking  in  the  polar  diagrams 
of  Figures  13  and  14,  derived  from  the  same  observations  as  Figure  12. 
Here  the  azimuthal  co-ordinate  is  System  III  longitude,  while  the 
radial  co-ordinate  is  the  number  of  intervals  of  millisecond  pulse 
activity  received  when  each  longitude  was  on  the  central  meridian. 
Figure  13  includes  only  those  intervals  when  the  System  III  longitude 
of  Io  lay  between  the  left-hand  pair  of  dashed  lines  in  Figure  12, 
that  is,  along  the  arc  shown  by  the  arrow  in  Figure  13.   Similarly, 
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FLA-CHILE    1963 
10  M/SLrC 


180 


NORTH    POLE 


270 


Figure  13.   Polar  plot  of  18  MHz  millisecond  pulse  activity  for  the 
indicated  range  of  Io  positions  centered  on  210°. 
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180 


FLA-CHILE      1963 
18  M/SEC 


NORTH    POLE 


270 


Figure  14.   Polar  plot  of  18  MHz  millisecond  pulse  activity  for  the 
indicated  range  of  Io  positions  centered  on  270°. 
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Figure  14  includes  only  data  recorded  when  Io  was  between  the  right- 
hand  pair  of  dashed  lines  in  Figure  12.   Figure  13  shows  that  when  B1 
is  active  on  the  central  meridian  (that  is,  directed  toward  the  earth), 
Io  is  near  System  III  longitude  210°,  and  also  roughly  90°  from  superior 
conjunction.   Figure  14  shows  that  millisecond  emission  from  B  also  is 
received  when  Io  is  about  90°  from  superior  conjunction,  but  in  this 
case  the  System  III  longitude  of  the  satellite  must  be  near  270°. 

Because  the  System  III  longitude  of  Io  advances  27.5°/hour, 
2.2  hours  are  required  for  Io  to  move  from  210°  to  270°,  but  during 
this  interval  the  longitude  of  the  central  meridian  increases  80°. 
Suppose  B,  is  on  the  central  meridian,  with  Io  at  System  III  longitude 
210°.   By  the  time  Io  has  moved  into  position  to  stimulate  B  ,  that 
source  is  already  well  past  the  central  meridian.   Although  the  peaks 
have  finite  widths,  Figures  13  and  14  indicate  that  the  overlap  is 
so  limited  that  it  is  improbable  that  both  B,  and  B  will  be  detected 
during  the  same  Jovian  rotation.   This  hypothesis  was  checked  by  an 
inspection  of  the  data  outputs  for  the  entire  apparitions  of  1963 
and  1964,  which  showed  that  on  no  occasion  were  millisecond  pulses 
received  from  both  B1  and  B  during  the  same  Jovian  rotation.   The 
same  was  generally  true  of  character  2  pulses,  although  on  several 
occasions  during  the  two  years  this  type  of  pulse  was  received  from 
B  and  B  during  the  same  rotation;  it  is  possible  that  these  exceptions 
represent  what  has  been  termed  "non-Io -related"  emission.   The  fact 
that  B  and  B  are  stimulated  by  distinctly  different  Io  positions 
is  perhaps  the  strongest  evidence  to  date  of  the  reality  of  the  division 
of  B  into  subsources . 

The  correlations  demonstrated  in  these  studies  can  be  used  to 
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support  the  hypothesis  that  the  immediate  environment  of  the  planet 
itself  must  contribute  in  an  important  way  to  the  burst  structure. 
This  point  of  view  received  additional  encouragement  from  Riihimaa's 
(1966)  report  that  the  bursts  from  source  B  display  a  distinctive 
spectral  fine  structure.   However,  the  foregoing  conclusion  should  not 
be  regarded  as  precluding  interplanetary  or  ionospheric  influences. 
For  example,  the  mechanism  responsible  for  the  generation  and  escape 
of  the  radiation  may  produce  among  the  several  sources  disparate 
angular  sizes  or  bandwidths .   These  factors,  then,  might  well  result 
in  the  radiation  from  each  source  interacting  differently  with  the 
intervening  medium. 


CHAPTER  IV 

VARIATIONS  OF  ROTATION  PERIOD  BY  SOURCE  AND  BY  FREQUENCY 

Introduction 

A  great  deal  of  interest  was  aroused  in  1962  when  a  change  in 
the  positions  of  the  Jovian  radio  sources  was  discovered  (Douglas 
and  Smith  1963b;  Smith  et_  al_.  1965a).   Since  most  of  the  theories 
of  the  decametric  radiation  linked  the  sources  to  the  Jovian  magnetic 
field,  which  in  turn  was  presumed  to  be  associated  with  the  planet's 
core,  this  drift  of  the  sources  implied  a  possible  change  in  the  rotation 
period  of  the  core. 

The  data  from  the  first  few  years  after  1960  suggested  that 
the  rotation  period  of  the  sources  had  suddenly  changed  from  the 
System  III  (1957.0)  period  to  one  1.17  seconds  longer  (Smith  et  al. 
1965a).   This  change  in  period  is  equivalent  to  a  drift  rate  of  10.45° 
per  year.   Lebo  (1964)  first  studied  this  drift  phenomenon  in  detail 
through  the  1962  apparition.   His  work  was  extended  to  include  the  1963 
apparition  and  greatly  expanded  in  scope  by  Hayward  (1965).   Upon 
Hayward's  departure  from  the  Florida  group,  the  study  of  the  drift  was 
continued  by  the  author  to  include  the  next  three  Jovian  apparitions 
(Olsson  and  Smith  1968) .   The  subject  of  this  chapter  is  the  study  of  the 
source  drift  through  the  1967.0  apparition.   More  recent  studies  that 
have  been  made  by  Register  (1968),  Donivan  and  Carr  (1969),  and  Carr 
et  al.  (1970)  include  one  more  year's  data. 
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Method  of  Analysis 

Histograms  of  the  reception  probability  versus  the  System  III 
longitude  of  the  central  meridian  were  generated  from  all  of  the  Florida 
and  Chile  data  at  the  frequencies  of  15,  18,  22.2,  and  27.6  MHz. 
The  Florida  and  Maipu,  Chile,  data  were  merged  into  a  single  histogram 
for  all  years  that  data  were  available  from  both  observatories.   Only 
these  merged  histograms  appear  in  the  appendix,  as  the  single -station 
histograms  have  been  reproduced  in  several  publications  and  dissertations 
by  the  Florida  group. 

The  method  used  to  determine  the  source  position  has  become 
somewhat  standardized  by  our  group.   A  tracing  paper  overlay  is  held 
in  place  on  the  histogram  and  one  approximates  the  source  as  well  as 
possible  with  a  roughly  Gaussian  curve.   The  longitude  of  the  peak  of 
this  sketched  curve  is  taken  as  the  location  of  the  source.   One  can 
easily  estimate  the  error  involved  by  slowly  shifting  the  curve  on  the 
tracing  paper  back  and  forth  in  longitude  until  the  fit  to  the  histogram 
is  obviously  poor.  At  frequencies  above  15  MHz  one  can  usually  locate 
the  sources  to  within  ±5°  using  this  method.  Whenever  the  source  was 
too  poorly  defined  or  too  asymmetric  to  use  the  above  method  accurately, 
the  source  center  was  defined  as  the  "centroid"  of  probability;  that 
is,  the  reception  probability  for  each  5°  zone  was  summed  over  the  entire 
source  and  the  longitude  of  the  midpoint  of  this  sum  was  determined. 
The  source  positions  are  listed  in  Table  3. 

Frequency  and  Source  Dependence  of  the  Drift 

Most  of  the  early  studies  concerning  the  drift  used  only  the 
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TABLE   3 
LOCATION  OF  THE   SOURCES 


Frequency 

Mean  Epoch 

Station* 

Source  B 

Source  A 

Source  C 

15  MHz 

1961.5 

M 

142± 

5° 

242± 

5° 

317±   5° 

1962.7 

F 

i20±: 

L0° 

260+: 

L0° 

325±   5° 

1963.8 

F 

150±: 

L0° 

255±: 

L0° 

325±   5° 

1964.9 

F 

140±: 

L0° 

250±: 

L0° 

315±   5° 

1965.9 

F 

165±: 

L0° 

260±! 

L0° 

328±   5° 

1967.0 

F 

166±: 

L0° 

260+: 

L0° 

330±   5° 

18  MHz 

1957.1 

F 

110± 

5° 

217± 

5° 

306±   5° 

1958.1 

F 

- 

23  5± 

5° 

309±   5° 

1959.1 

F 

- 

215± 

5° 

__ 

1960.2 

FM 

125± 

5° 

230± 

5° 

323±   5° 

1961.6 

FM 

130± 

5° 

23  5± 

5° 

312±   5° 

1962.7 

FM 

130+ 

5° 

249  ± 

5° 

320±   5° 

1963.8 

FM 

141± 

5° 

260± 

5° 

327±   5° 

1964.8 

FM 

137± 

5° 

268± 

5° 

334±   5° 

1965.9 

FM 

156± 

5° 

265± 

5° 

325±   5° 

1967.0 

FM 

160± 

5° 

260± 

5° 

332±   5° 

22.2  MHz 

1958.1 

F 

105±: 

L0° 

220± 

5° 

__ 

1960.3 

FM 

ii4±: 

L0° 

220± 

5° 

315±   5° 

1961.5 

FM 

132± 

5° 

227± 

5° 

299±   5° 

1962.8 

F 

130± 

5° 

243± 

5° 

312±   5° 

1963.8 

FM 

140± 

5° 

255± 

5° 

320±   5° 

1964.8 

FM 

138± 

5° 

257± 

5° 

320±   5° 

1965.9 

FM 

143± 

5° 

247± 

5° 

315±   5° 

1967.0 

FM 

147± 

5° 

240± 

5° 

319±   5° 

27.6  MHz 

1961.6 

F 

115± 

5° 

215± 

5° 

298±   5° 

1962.7 

F 

125± 

5° 

240± 

5° 

335±   5° 

1963.7 

FM 

133± 

5° 

233± 

5° 

327±   5° 

1964.7 

FM 

140± 

5° 

245± 

5° 

297±   5° 

1965.9 

F 

140± 

5° 

245± 

5° 

318±   5° 

1967.0 

F 

148± 

5° 

232± 

5° 

~- 

*  In  the  station  column,  F   =  Florida   and  M  =  Maipu 
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position  of  source  A  as  determined  from  the  data  collected  at  frequencies 
from  18  to  20  MHz,  the  reason  being  that  A  is  the  most  prominent  source 
and  the  data  at  these  frequencies  are  the  most  abundant.   The  question 
of  frequency  and  source  dependence  of  the  drift  was  investigated  by 
Hayward  (1965).   He  found  that  the  drift  rate  averaged  over  a  three  - 
year  period  did  vary  with  frequency  in  a  fairly  smooth  fashion  -- 
the  higher  the  frequency,  the  lower  the  drift  rate.   This  same  trend 
was  still  apparent  with  the  inclusion  of  another  year  of  data  (Smith 
et  al.  1965b).   Analysis  of  the  next  two  years'  data  showed  that  some 
of  the  sources  had  reversed  their  sense  of  drift  and  were  then  migrating 
toward  lower  longitudes.   A  study  of  all  of  the  data  by  Gulkis  and 
Carr  (1966)  revealed  that  source  A  may  be  undergoing  a  sinusoidal 
drift  in  longitude  with  a  period  of  about  12  years  . 

To  investigate  the  year -by -year  drift  of  the  three  sources  A, 
B,  and  C  at  the  frequencies  of  15,  18,  22.2,  and  27.6  MHz,  the  longitude 
of  each  source  was  plotted  as  a  function  of  date  for  each  frequency. 
Figures  15  through  17  resulted.   In  each  of  the  figures  the  18  MHz 
data  are  included  as  a  comparison.   Error  bars  were  not  shown  for  those 
pairs  of  data  points  where  the  separation  was  less  than  the  error.   Lines 
were  drawn  joining  adjacent  points  of  a  given  frequency,  and  the  slopes, 
which  are  the  average  yearly  drift  rates,  were  calculated  and  recorded 
in  Table  4. 

Figure  15  shows  that  all  three  of  the  sources  at  15  and  18  MHz 
drifted  together  within  the  limits  of  the  accuracy  of  the  measurement 
for  all  years,  except  for  sources  A  and  B  in  1964.8.   Figure  16  shows 
that  the  three  sources  at  18  and  22.2  MHz  drifted  together  until  1964.8, 
when  sources  A  and  C  began  to  separate,  followed  by  source  B  a  year  later. 
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TABLE  4 
YEARLY  AVERAGE  DRIFT  RATE  OF  SOURCES 
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Frequency- 

Mean 

Station 

Source   B 

Source  A 

Source  C 

Epoch 
1962.7 

Drift  Rate 
-18.3 

Drift  Rate 
15.0 

Drift  Rate 

lS  MHz 

F 

6.7 

1963.8 

F 

27.3 

-4.6 

.0 

1964.9 

F 

-9.1 

-4.6 

-9.1 

1965.9 

F 

25.0 

+10.0 

+13.0 

1967.0 

F 

0.9 

.0 

1.8 

18  MHz 

1958.1 

F 

__ 

18.0 

3.0 

19  59.1 

F 

-- 

-20.0 

-- 

1960.2 

FM 

4.8 

13.6 

6.7 

1961.6 

FM 

4.6 

3.6 

-7.9 

1962.7 

FM 

.0 

12.7 

7.3 

1963.8 

FM 

10.0 

10.0 

6.4 

1964.8 

FM 

-4.0 

8.0 

7.0 

1965.9 

FM 

17.3 

-2.7 

-8.2 

1967.0 

FM 

3.6 

-4.6 

6.4 

22.2  MHz 

1960.3 

FM 

4.1 

.0 

__ 

1961.5 

FM 

-16.7 

5.8 

-13.3 

1962.8 

F 

-1.5 

12.3 

10.0 

1963.8 

FM 

10.0 

12.0 

8.0 

1964.8 

FM 

-2.0 

2.0 

.0 

1965.9 

FM 

4.5 

-9.1 

-4.5 

1967.0 

FM 

3.6 

-6.4 

3.6 

27.6  MHz 

1962.7 

F 

9.1 

22.7 

33.6 

1963.7 

FM 

7.3 

-6.4 

-7.3 

1964.8 

F 

6.4 

10.9 

-27.3 

1965.9 

F 

.0 

.0 

19.1 

1967.0 

F 

7.3 

-11.8 

— 
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The  most  striking  feature  in  Figure  17  is  the  wide  separation  between 
source  A  at  18  and  27.6  MHz  for  all  apparitions  except  1962.7.   In 
these  figures  there  are  57  pairs  of  points.   For  33  of  them  the  sepa- 
ration is  less  than  the  error  bars,  and  for  the  remaining  24  the 
separation  is  greater.   It  should  be  noted  that  the  separations  were 
a  minimum  in  all  of  the  figures  for  the  three -year  period  from  1962.7 
to  1964.8,  when  the  amount  of  data  was  far  greater  than  for  apparitions 
before  or  after.   This  suggests  that  the  deviations  observed  for  other 
years  may  be  the  result  of  poorer  statistics. 

Figures  18  through  22  were  drawn  to  illustrate  the  differences 
in  the  drift  rates  from  source  to  source  at  a  fixed  frequency  and  also 
the  differences  in  the  drift  rates  for  a  given  source  at  different 
frequencies.   Figures  18  and  19  are  plots  of  the  drift  rate  of  the 
three  sources  as  a  function  of  frequency  for  the  apparitions  of  1964.8 
and  1965.9,  respectively.  Note  in  these  two  figures  that  a  point 
below  the  zero -axis  represents  a  negative  drift  (that  is,  the  source 
moved  to  a  lower  longitude  during  the  year) ,  and  a  point  above  the 
zero -axis  represents  a  positive  drift.  These  two  figures  show  that 
the  drift  rate  is  in  the  same  sense  for  all  three  sources  only  at 
15  MHz.   For  all  of  the  other  frequencies  the  sense  of  drift  is  different 
for  at  least  two  of  the  sources . 

Figures  20  through  22  are  plots  of  the  drift  rate  of  the  sources 
as  a  function  of  time  for  the  four  different  frequencies.  These  figures 
show  that  there  is  a  wide  variation  in  the  magnitude  and  direction  of 
the  drift  rate  of  the  sources  for  the  different  frequencies  over  the 
period  of  time  indicated.   There  seems  to  be  little,  if  any,  correlation 
between  the  frequency  and  the  drift  rate  of  a  particular  source  for  a 
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Figure  18.   Source  drift  rate  versus  frequency  for  1964.9, 
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Figure  19.   Source  drift  rate  versus  frequency  for  1963.9. 
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Figure  20.   Drift  rate  of  source  B  as  a  function  of  time, 
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Figure  21.   Drift  rate  of  source  A  as  a  function  of  time, 
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Figure  22.   Drift  rate  of  source  C  as  a  function  of  time, 
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given  year.   The  data  in  these  figures  suggest  that  the  distance  between 
the  sources  would  also  vary  year  by  year.   They  also  suggest  a  yearly 
variation  in  the  distance  between  the  positions  of  any  given  source 
measured  at  two  different  frequencies.   To  check  for  this  variation, 
the  difference  in  position  of  source  A  as  determined  at  18  and  15  MHz 
was  calculated  for  each  year.   These  differences  were  also  calculated 
for  source  A  using  the  frequencies  18  and  22.2  MHz  and  18  and  27.6  MHz. 
Similar  differences  were  calculated  for  sources  B  and  C.   These  results 
are  given  Table  5.  All  of  the  differences  are  absolute  values. 
An  average  difference  was  calculated  for  each  column  in  the  table  and  the 
error  was  computed.  The  error  was  assumed  to  be  equal  to  the  square 
root  of  the  sum  of  the  squares  of  the  error  of  each  source  position. 
Differences  between  two  adjacent  sources  at  a  fixed  frequency  were 
also  calculated  for  each  year  and  are  listed  in  Table  6. 

Of  the  57  longitude  differences  recorded  in  Table  5,  only  5 
(or  about  9%)  fall  outside  of  the  computed  error.   This  suggests 
that  the  year -by -year  source  location  is  independent  of  frequency. 
It  should  be  understood  that  this  result  is  not  in  conflict  with  an 
earlier  study  by  the  Florida  group  of  a  frequency -dependent  drift  of 
the  sources  (Smith  et_  al_.  1965a)  .   In  the  earlier  study,  the  drift 
was  the  change  in  the  locations  of  the  sources  with  a  change  in  frequency 
for  a  given  apparition,  as  shown  in  Figure  23,  which  is  copied  from 
the  paper  by  Smith  et_  al_.   The  present  work  investigated  the  differences 
in  source  position,  such  as  the  interval  indicated  in  Figure  23  by 
the  line  element  labeled  d,  and  noted  how  these  differences  change 
from  one  apparition  to  the  next. 

In  Table  6,  21  of  the  54  differences  (or  about  39%)  fall  outside 
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TABLE  5 

LONGITUDE  DIFFERENCES  IN  SOURCE  LOCATION 
AT  DIFFERENT  FREQUENCIES 


Mean  Epoch 


19  58 
1960 


A18-A22* 

15° 
10° 


A15-A18 


A18-A27 


1961.5 

8° 

7° 

20° 

1962.8 

6° 

11° 

9°** 

1963.8 

5° 

5° 

27° 

1964.8 

11° 

18° 

23° 

1965.9 

18° 

5° 

20° 

1967.0 

20° 
Av.=   12±7° 

Av 

0° 
.=  8±14° 

Av 

28° 
.=  21±7° 

B18-B22 

] 

315-B18 

: 

B27-B18 

1960.3 

11° 

__ 

_. 

1961.5 

2° 

12° 

15° 

1962.8 

0° 

10° 

5° 

1963.8 

1° 

9° 

8° 

1964.8 

1° 

3° 

3° 

1965.9 

13° 

9° 

16° 

1967.0 

13° 
Av.=  6±7° 

Av 

6° 
.=  8±14° 

Av 

12° 
.=  10±7° 

C18-C22 

C15-C18 

C27-C18 

1960.3 

8° 

__ 

.. 

1961.5 

13° 

5° 

14° 

1962.8 

8° 

5° 

15° 

1963.8 

7° 

2° 

0° 

1964.8 

14° 

19°** 

37°** 

1965.9 

10° 

3° 

7°** 

1967.0 

13° 

2° 

— 

Av.=   10±7° 

Av 

.=  6±7° 

Av 

.=   15±7° 

*  The  notation  A18  represents  the  location  of  source  A  measured 
at  18  MHz,  etc. 
**  These  values  exceed  the  average  by  an  amount  greater  than 
the  error  . 


TABLE  6 

DIFFERENCES  IN  LONGITUDES  OF  THE  SOURCES 
AT  FIXED  FREQUENCY 
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Frequency 


Mean  Epoch 


|C-A| 


|A-B| 


15  MHz 

1961.5 

75° 

100° 

1962.7 

65° 

140°* 

1963.8 

70° 

105° 

1964.9 

65° 

110° 

1965.9 

68° 

95° 

1967.0 

Av. 

70° 

Av. 

94° 

=  69±11° 

=   107±14° 

18  MHz 

1957.1 

89 1* 

— 

1958.1 

74° 

-- 

1960.2 

93°* 

105°* 

1961.6 

77° 

105°* 

1962.7 

71° 

119° 

1963.8 

67° 

119° 

1964.8 

66°* 

131°* 

1965.9 

60°* 

109° 

1967.0 

Av. 

72° 

Av. 

100°* 

=  74±7° 

=   113±7° 

22.2  MHz 

19  58.1 

__ 

115°* 

1960.3 

95°* 

106° 

1961.5 

72° 

95°* 

1962.8 

69° 

113° 

1963.8 

65°* 

115°* 

1964.8 

63°* 

119°* 

1965.9 

68° 

104° 

1967.0 

Av 

79° 

Av 

93°* 

=  73±7° 

=   107±7° 

27.6  MHz 

1961.6 

83° 

100° 

1962.7 

95°* 

115° 

1963.7 

94°* 

100° 

1964.7 

52°* 

105° 

1965.9 

73° 

105° 

1967.0 

Av 

-- 

Av 

84°* 

=  79±7° 

=   101±7° 

*These  values  exceed  the  average  by  an  amount  greater  than 
the  error  . 


Figure  23.   Positions  of  sources  A,  B,  and  C  at  various 
frequencies  in  1962,  with  lines  tracing  a 
possible  mode  of  "drift"  of  the  sources. 
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of  the  calculated  error,  indicating  a  possible  variation  in  the  distances 
between  the  sources.   This  variation  could,  however,  be  due  to  the 
difficulty  one  has  in  accurately  locating  sources  B  and  C  for  years 
of  low  activity. 


CHAPTER  V 

A  STUDY  OF  THE  ROTATION  PERIOD  BASED 
ON  STORM  COMMENCEMENT  TIMES 

Introduction 


Several  hypotheses  have  been  suggested  to  explain  the  apparent 
sinusoidal  drift  of  source  A.   Gulkis  and  Carr  (1966)  assume  that  the 
radiation  from  the  source  is  beamed  into  a  thin,  curved  sheet  which 
is  fixed  with  respect  to  the  Jovian  magnetic  field.   Since  Jupiter's 
magnetic  axis  is  known  to  be  tilted  some  10°  with  respect  to  its 
rotational  axis,  the  sheet  of  radiation  will  wobble  as  viewed  from  the 
earth  as  the  planet  rotates .   Thus  one  can  receive  radiation  only  for 
certain  orientations  of  the  earth  and  Jupiter's  magnetic  field.   This 
hypothesis  therefore  explains  the  variation  in  source  position  in 
terms  of  the  rate  of  change  of  the  Jovicentric  declination  of  the 
earth.   The  hypothesis  has  been  strengthened  by  the  recent  work  of 
Donivan  and  Carr  (1969)  and  Carr  et_  al_.  (1970),  who  have  included 
several  more  years  of  data  in  the  study,  which  continues  to  show  a 
strong  correlation  between  source  position  and  the  Jovicentric  decli- 
nation of  the  earth. 

Another  possible  explanation  of  the  drift  of  source  A  was 
recently  proposed  by  Register  and  Smith  (1969).   For  several  years 
it  has  been  known  that  source  A  seems  to  be  composed  of  an  Io -related 
component  and  a  non-Io -related  component  (Olsson  and  Smith  1966; 
Shever  1967)  .   If  one  separates  the  radiation  from  the  source  into 
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these  two  components  and  plots  yearly  histograms  from  these  data, 
several  differences  are  noted  between  the  curves.   The  centers  of 
source  A  of  the  Io -related  curves  consistently  occur  at  a  lower 
longitude  than  do  the  non-Io-related  centers.   Secondly,  the  source 
centers  display  completely  different  drift  characteristics.   The  drift 
of  the  center  of  the  source  in  the  Io -related  curves  is  almost  linear, 
whereas  the  non-Io-related  curves  strongly  suggest  a  sinusoidal  drift. 
These  histograms  also  show  that  while  the  amplitude  of  the  Io-related 
source  is  fairly  constant  year  after  year,  the  non-Io-related  source 
varies  considerably  in  amplitude.   Therefore,  it  is  possible  that 
the  drift  of  the  blended  source  is  the  result  of  superposing  a  non-Io- 
related  component  of  variable  amplitude  and  variable  period  on  a  more 
stable  Io-related  component. 

Duncan's  Hypothesis 

A  third  explanation  has  been  offered  by  Duncan  (1967) .   Approaching 
the  question  from  a  different  point  of  view,  Duncan  defined  the  Jovian 
radio  rotation  period  in  terms  of  the  commencement  times  of  the  storms 
from  source  A.   Duncan's  analysis  utilized  the  observations  of  other 
workers,  which  contained  just  over  1000  storms  from  all  Jovian  sources 
at  frequencies  ranging  from  18  MHz  to  30  MHz.   He  plotted  histograms 
of  storm  commencement  times  and  showed  that  the  leading  edge  or  front 
of  source  A  seemed  to  remain  stationary  in  a  coordinate  system  rotating 
with  a  period  of  9  55  29.70.  When  he  studied  the  center  of  the  source 
in  these  histograms,  Duncan  found  that  some  source  drift  remained, 
although  he  regarded  it  as  considerably  smaller  than  that  reported  by 
others.   He  concluded  from  these  observations  that  the  source  drift 
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was  in  fact  due  entirely  to  a  long-term  secular  variation  in  the 
average  durations  of  the  Jovian  noise  storms ,  an  effect  that  would 
shift  the  center  of  the  apparent  source  but  not  its  leading  edge. 

Analysis  of  Florida  and  Chile  Observations 

Since  the  extensive  data  of  the  University  of  Florida  and 
University  of  Chile  radio  observatories  were  not  available  to  Duncan, 
it  seemed  worth  while  to  examine  them  as  a  test  of  his  very  interesting 
hypothesis  (Olsson  et  al .  1969a,  1969b).   The  study  was  confined  to 
the  18  and  22.2  MHz  observations,  which  are  consistently  the  most 
numerous  and  most  reliable.   The  University  of  Florida  observations 
utilized  in  this  investigation  extend  over  ten  apparitions,  while  those 
from  Chile  cover  seven  apparitions.   The  daily  records  from  each 
station  were  examined  and  the  System  III  (1957.0)  central  meridian 
longitude  was  recorded  for  the  time  at  which  each  noise  storm  began. 
Groups  of  bursts  separated  by  intervals  of  quiet  of  less  than  15 
minutes  (about  9°  of  Jovian  rotation)  were  regarded  as  belonging 
to  the  same  storm.  Altogether,  over  1800  noise  storms  from  source  A 
were  included  in  the  analysis.   A  histogram  was  then  prepared  for 
each  apparition,  showing  the  number  of  storm  commencements  that  had 
been  observed  when  each  5°  zone  of  System  III  longitude  was  crossing 
the  central  meridian.   Since  it  is  known  that  the  source  structure  in 
the  ordinary  probability  histograms  changes  with  frequency  (Smith 
et_  al_.  1965a)  ,  separate  histograms  were  plotted  for  18  MHz  and  for 
22.2  MHz.   To  provide  a  cross-check,  individual  histograms  were  also 
prepared  for  each  of  the  two  observatories. 

The  33  histograms  that  resulted  from  this  analysis  are  shown  in 
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Figures  24  through  28.   It  should  be  noted  in  several  of  the  histograms 
that  the  source  A  peak  is  bifurcated;  that  is,  there  is  a  second 
grouping  of  commencements  at  the  higher  longitudes .   Those  histograms 
which  display  this  bifurcation  most  prominently  are  the  18  MHz  Florida 
ones  for  1963.8  to  1965.9,  the  22  MHz  Florida  ones  for  1964.8  to  1967.1, 
and  the  22  MHz  Chile  histogram  for  1965.9.   This  same  bifurcation  is 
evident  in  Duncan's  second  histogram  (Duncan  1967,  Figure  6),  which 
is  a  merge  of  data  spanning  the  years  1962-1966,  although  it  is  absent 
in  his  earlier  histogram  for  the  period  1951-1961.   Figure  29  is  a 
merge  of  the  18  MHz  Florida  data  for  the  periods  1958-1961  and  1962- 
1966.   These  data  produce  a  distribution  strikingly  similar  to  Duncan's 
histograms.   In  both  cases  the  steepest  part  of  the  commencement 
front  is  preceded  by  a  more  gradual  rise  covering  about  25°  in  longitude 
The  bifurcated  peaks  span  approximately  60°  and  are  followed  by 
tails  extending  to  longitudes  near  300°. 

Figures  30  through  33  are  the  histograms  produced  when  the 
Florida  and  Chile  data  were  merged.   For  each  curve  the  "front"  was 
selected  as  the  first  significant  rise  following  the  well-defined 
minimum  which  characteristically  separates  sources  A  and  B.   It  is 
evident  from  the  method  of  analysis  that  any  movement  of  the  front 
must  occur  in  multiples  of  5°.   Furthermore,  in  several  instances 
there  was  more  than  one  feature  that  might  have  been  selected  as  the 
front;  where  these  ambiguities  exist  two  choices  were  made.   The 
fronts  so  chosen  are  indicated  on  the  histograms  of  the  merged  data 
with  vertical  dashed  lines . 

Figure  34  traces  the  position  of  the  commencement  front  of  source 
A  as  it  was  determined  from  the  merged  Florida -Chile  histograms. 
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Figure  24.   Commencement  time  of  storms  histograms  of 
source  A  for  1965.9  to  1968.1  at  18  MHz. 
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Figure  25.   Commencement  time  of  storms  histograms  of 
source  A  for  1958.1  to  1964.8  at  18  MHz. 
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Figure  26.   Commencement  time  of  storms  histograms  of 
source  A  for  1960.2  to  1967.0  at  18  MHz. 


66 


FLORIDA   2?. 2  MHz SOURCE    A 

'              I  '  h 

1968.!        (-LJ L_ T 


0 
10 


Lr 


£ JLL 


19  STORMS 
41  STORMS 


150  200"  2150"  300"  350v 

LONGITUDE    OF   CENTRAL    MERIDIAN 


Figure  27.   Commencement  time  of  storms  histograms  of 
source  A  for  1958.1  to  1968.1  at  22.2  MHz, 
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Figure  28.   Commencement  time  of  storms  histograms  of 
source  A  for  1960.2  to  1967.0  at  22.2  MHz. 
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Figure  29.   Merged  histograms  of  commencement  time  of  storms 
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Figure  30.   Histograms  showing  commencement  fronts  with 
dashed  lines  for  1964.8  to  1968.1  at  18  MHz. 
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Figure  31.   Histograms  shoving  commencement  fronts  with 
dashed  lines  for  1958.1  to  1963.8  at  18  MHz, 
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Figure  32.   Histograms  showing  commencement  fronts  with 

dashed  lines  for  1964.8  to  1968.1  at  22.2  MHz, 
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Figure  33.   Histograms  showing  commencement  fronts  with 

dashed  lines  for  19  58.1  to  1963.8  at  22.2  MHz. 
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Figure  34.   Location  of  the  commencement  front  of  source  A 
during  ten  apparitions  of  Jupiter. 
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Where  two  fronts  were  chosen  for  one  year,  both  are  recorded  and  joined 
by  a  horizontal  tie.   The  straightforward  interpretation  of  Figure  34 
is  that  it  supports  Duncan's  contention;  over  an  interval  of  a  decade 
the  commencement  front  has  remained  substantially  fixed  in  longitude, 
which  is  to  say  that  its  rotation  period  has  been  constant.   It  is 
also  interesting  to  note  that  the  longitude  of  the  front  is  the  same 
at  18  and  22.2  MHz,  whereas  there  is  generally  a  significant  displacement 
of  the  center  of  the  source  between  these  frequencies  in  the  conventional 
probability  histograms. 

The  fact  that  each  set  of  points  in  Figure  34  is  well  fitted 
by  a  vertical  straight  line  implies  that  the  period  defined  by  the 
commencement  fronts  in  this  study  is  that  of  System  III  (1957.0), 

i.e.,  9  55m29f37  with  a  standard  error  of  ±  o!o8.   This  is  about 

s 
0.3  less  than  the  period  favored  by  recent  decimetric  observations, 

as  well  as  the  mean  period  adopted  by  Gulkis  and  Carr  (1966)  in  their 

oscillatory  theory  of  the  behavior  of  the  standard  probability  histograms. 

In  the  present  state  of  the  art  it  is  well  to  exercise  caution  in 

assuming  that  each  of  these  highly  disparate  approaches  is  measuring 

precisely  the  same  phenomenon. 

Figure  35  shows  the  location  of  the  center  of  source  A  in  the 

storm  commencement  histograms  for  the  ten  apparitions.   The  points 

were  obtained  by  determining  the  central  meridian  longitude  corresponding 

to  the  average  storm  commencement  time  for  each  apparition.   It  is 

evident  that,  unlike  the  commencement  front,  the  source  center  displays 

an  oscillatory  drift  of  large  amplitude;  the  solid  curves  are  sinusoids 

with  11. 9 -year  periods  which  were  least -squares  fitted  to  the  empirical 

points.   The  dashed  curve,  adapted  from  Gulkis  and  Carr  (1966),  depicts 
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the  oscillation  of  the  center  of  source  A  in  the  usual  probability 
histograms.   It  is  evident  that  the  drift  amplitude  of  the  source 
center  is  substantially  the  same  in  both  types  of  histograms. 

Storm  Durations 


A  final  question  is  whether  the  oscillatory  behavior  in  Figure  3  5 
can  be  attributed  to  a  long-term  periodicity  in  the  mean  durations 
of  the  noise  storms .   Experienced  observers  are  aware  that  at  any 
frequency  it  is  unusual  for  a  storm  simply  to  begin  at  the  leading 
edge  of  the  source  (point  a_  in  Figure  36)  and  persist  right  across 
the  source  to  b_.  Far  more  commonly,  storms  begin  "within"  the  source 
at  some  point  such  as  c_  and  end  at  another  point  d_  well  before  the 
trailing  edge  of  the  source  reaches  the  central  meridian.   In  other 
words,  the  "source"  is  largely  a  superposition  of  numerous  fairly 
short  storms  with  staggered  commencements  and  endings  (such  staggering 
of  at  least  the  commencements  is  of  course  implicit  in  the  distribution 
shown  in  Figures  24  through  28) . 

To  investigate  the  distribution  of  storms  within  source  A,  plots 
in  which  each  storm  is  represented  by  a  horizontal  line  were  prepared 
for  each  apparition  from  1959.1  to  1967.0  from  the  18  MHz  and  22.2  MHz 
Florida  data.  These  lines  are  plotted  against  the  longitude  of  the 
central  meridian,  and  thus  indicate  where  the  storm  began  and  how  long 
it  lasted.  The  results  are  shown  in  Figures  37  through  46.   In  each 
case  there  was  a  wide  distribution  of  storm  durations  and  of  starting 
and  stopping  times.   To  study  this  result  further,  a  distribution  curve 
was  plotted  as  a  function  of  the  duration  of  the  storms  and  this  is 
shown  in  Figure  47.   In  this  figure  the  duration  is  measured  in  5°  zones 
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Figure  37.   Distribution  of  storms  within  source  A  for 
1959.1  to  1961.6  at  18  MHz. 


80 


200  250"  300"  350^ 

LONGITUDE    OF  TME    CENTRAL    MERIDIAN 


Figure  38.   Distribution  of  storms  within  source  A  for 
1962.7  at  18  MHz. 
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Figure  39 .  Distribution  of  storms  within  source  A  for 
1963.8  at  18  MHz. 


82 


150      200      250      300       350 
LONGITUDE  OF  THE  CENTRAL  MERIDIAN 


Figure  40.   Distribution  of  storms  within  source  A  for 
1964.8  at  18  MHz. 
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Figure  41.   Distribution  of  storms  within  source  A  for 
1965.9  at  18  MHz. 
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Figure  42. 


Distribution  of  storms  within  source  A  for 
1967.0  at  18  MHz. 
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Figure  43.   Distribution  of  storms  within  source  A  for 
1959.1  to  1962.7  at  22.2  MHz. 


Sb 


2  2  MHz  FLA 

1963.8 
102    STORMS 


I 
150°  200°  250°  300°  350' 

LONGITUDE   OF    THE    CENTRAL    MERIDIAN 


Figure  44.   Distribution  of  storms  within  source  A  for 
1963.8  at  22.2  MHz. 
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Figure  45.   Distribution  of  storms  within  source  A  for 
1964.8  at  22.2  MHz. 
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Figure  46.   Distribution  of  storms  within  source  A  for 
1965.9  to  1967.0  at  22.2  MHz. 
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TABLE   7 

STORM  STATISTICS 

Apparition 

Frequence 

Number   of   Storms 

Average  Duration 

(MHz) 
18 

(Minutes) 

1959.1 

18 

37.2 

22 

6 

34.2 

1960.2 

18 

23 

36.6 

22 

10 

39.0 

1961.6 

18 

100 

73.8 

22 

31 

39.0 

1962.7 

18 

129 

52.8 

22 

51 

48.0 

1963.8 

18 

149 

78.0 

22 

102 

64.8 

1964.8 

18 

143 

52.2 

22 

84 

47.4 

1965.9 

18 

142 

47.4 

22 

81 

35.4 

1967.0 

18 

88 

45.0 

22 

53 

39.6 

91 


1968.0 


18  MHz 


22.2  MHz 


1966.0 


Id 
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Figure  48.      A  history  of  mean  storm  duration. 
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of  longitude,  where  1  zone  =8.27  minutes,  and  the  distribution  is 
expressed  as  the  fraction  of  the  total  number  of  storms  that  have  a 
particular  duration.   If  one  ignores  the  first  point  in  the  figure, 
the  others  are  best  fitted  by  a  monotonically  decreasing  curve  as 
shown.   There  is  justification  for  omitting  the  first  point,  which 
represents  storm  durations  from  a  millisecond  burst  up  to  storms 
lasting  8.27  minutes,  for  without  doubt  many  of  the  very  short  single 
bursts  have  either  been  missed  or  misidentif ied  by  the  observers; 
any  experienced  Jovian  observer  would  attest  to  this  statement. 
Therefore,  it  is  conceivable  that  the  first  point  in  the  figure 
represents  only  a  fraction  of  the  true  value.   Very  few  storms  with 
durations  greater  than  several  minutes  would  be  missed,  and  thus  the 
other  points  should  be  accurate. 

To  quantify  the  study,  the  mean  storm  duration  was  computed  for 
each  frequency  and  each  apparition,  with  the  results  shown  in  Table  7 
and  Figure  48.   At  neither  frequency  do  the  points  suggest  a  monotonic 
increase  in  storm  duration  from  1959.1  to  1967.0.   There  is  a  suggestion 
of  a  maximum  in  storm  duration  around  1963,  but  this  falls  mid -way 
between  the  positive  and  negative  peaks  of  the  sinusoids  in  Figure  35. 
Thus  the  analysis  does  not  appear  to  support  the  view  that  the  oscillation 
of  the  source  center  is  due  to  a  simple  cyclic  variation  in  storm 
duration. 


CHAPTER  VI 

SATELLITE  INFLUENCES  ON  THE  RADIATION 

Hypothetical  Satellites 

In  1967  Boswell  and  Wilson  (1967)  discovered  that  the  eight 
inner  satellites  of  Saturn  obey  the  following  rule: 

Ni+Ni+2=    <1+Mik>    (Ni+1+N.+3) 

where  N  is  the  mean  orbital  angular  rate  of  the  ith  satellite,  M. 
1  i 

is  a  small  positive  integer  and  k  is  a  constant.   Because  of  this 
relationship  they  reasoned  that  such  a  rule  certainly  does  not  arise 
by  chance.   They  then  found  two  solutions  of  the  equation  for  the  Jovian 
system  which  predicted  the  possible  existence  of  one  or  two  undiscovered 
satellites  orbiting  Jupiter  between  Io  and  Amalthea.   These  two  solutions 
gave  the  values  of  k  =  0.365514,  N  =  432.189°/day  and  k  =  0.262087, 
N  =  417 .802°/day .   No  values  were  given  for  M. .   The  sidereal  periods 
of  the  conjectured  satellites  were  given  as  20  40  27  and  19  59  24  . 
These  periods  correspond  to  mean  orbital  radii  of  3.66  and  3.58 
planetary  radii  respectively. 

This  prediction  motivated  Bigg  (1969)  to  investigate  the  decametric 
radiation  from  Jupiter  for  possible  effects  due  to  the  outermost 
satellite  of  Boswell  and  Wilson.   Bigg  used  the  Boulder  radio  spectrograph 
data  for  the  years  1964  to  1968  and  the  19.7  MHz  data  from  CSIRO  from 
1959  to  1963.   Bigg  reported  a  very  strong  modulation  of  the  emission 


93 


94 


at  a  period  of  20  40  57.4,  which  is  within  5  seconds  of  the  mean 
synodic  period  of  the  predicted  satellite.  He  also  found  a  modulation 
associated  with  a  period  of  20  43  36.8,  a  value  far  from  that  of  the 
second  predicted  satellite.   He  called  the  satellite  with  the  shorter 
period  "A"  and  the  other  "B". 

Reception  probability  was  plotted  versus  the  longitudes  of  A 
and  B,  and  diagrams  quite  similar  to  the  Io  curves  resulted.   The  influence 
of  A  was  similar  to  Io's  but  not  as  strong.   The  B  satellite  plot  had 
only  one  peak,  rather  than  two  as  in  the  case  of  Io .   These  results 
convinced  Bigg  that  there  is  a  strong  possibility  that  satellites  A 
and  B  exist.   One  should  note  that  Bigg's  satellite  B  is  not  the  same 
as  Boswell  and  Wilson's  second  one.  With  further  analysis,  he  was  able 
to  calculate  the  longitudes  of  the  two  satellites  at  0  hours  U.  T.  on 
December  31,  1959.   They  were:  A,  340°  and  B,  180°. 

The  Florida  Analysis 

Once  again  it  was  felt  that  Bigg's  discovery  of  a  satellite 
influence  on  the  Jovian  emission  should  be  checked  using  the  more 
extensive  Florida -Chile  data.   This  check  was  accomplished  very  easily 
by  using  the  newly  developed  computer  program  (Register  1968) .   An 
ephemeris  of  the  position  from  superior  conjunction  at  0  hours  U.  T. 
for  each  day  was  generated  for  A  and  B.   These  data  were  put  into  the 
new  program  and  a  probability  histogram  for  the  18  MHz  data  for  the  1963.8 
Florida -Maipu  apparition  was  obtained.   Figure  49  is  the  histogram, 
showing  the  probability  as  a  function  of  satellite  A's  position. 
Three-point  smoothing  was  used  in  plotting  this  curve  to  reduce  statistical 
fluctuations.  At  first  glance,  the  figure  suggests  two  broad  peaks 


95 


o 
O 
'O 

to 


o 

o 

CD 

rO 


O 

H 
O 
Z. 

=> 
"3 

o 
o 


°  o 
o  — 

c\j  q: 

i.j 

Q. 

ID 
CO 

o      ^ 

u)  O 
-_  CC 

u. 


Ui 
O  -I 

—  _J 

LU 

h- 

< 
en 


o 
O 


Ainiavaoad 


96 


separated  by  a  narrow  valley  at  about  180°.   However,  when  the  same 

data  were  replotted  with  180°  at  the  origin  of  the  horizontal  axis, 

it  appears  that  there  is  but  one  broad  peak.   Figure  50  is  the  plot 

with  the  180°  shift  of  the  coordinates.   Thus  the  most  prominent 

feature  is  the  decrease  in  probability  when  A  is  near  inferior  conjunction. 

This  feature  also  appeared  when  the  Florida  and  Maipu  data  were  plotted 

separately.   The  histogram  for  satellite  B  was  similar  to  Figure  49 

with  a  far  less  pronounced  decrease  in  reception  probability.  As  a 

further  check  for  possible  A  and  B  influences,  all  of  the  raw  data 

computer  cards  containing  Io -related  data  were  removed  from  the  deck 

and  the  program  was  rerun.   Figure  51  is  the  resulting  histogram. 

Even  with  the  Io -related  data  omitted,  there  is  no  clear-cut  relation 

between  the  probability  and  the  position  of  A. 

In  a  paper  which  appeared  after  the  above  analysis,  Douglas  and 
Bozyan  (1970)  presented  a  convincing  argument  refuting  Bigg's  work  by 
using  a  different  approach  than  that  described  above.   They  showed 
that  since  the  Jovian  radiation  and  the  detection  of  the  radiation 
depend  upon  at  least  four  periods,  harmonics  of  these  periods  as 
well  as  their  sums  and  differences  will  appear  in  the  data.   They 
showed  that  one  would  expect  spurious  periodicities  to  be  present 
with  mean  synodic  periods ,  P ,  given  by 

P"1  =  jPj1  +  kP'J  +  mP'1  +  nP~J  ,  j,  k,  m,  n  =  0,  ±1,  ±2,  ... 

where  P  ,  P   ,  P   are  the  mean  synodic  periods  of  Jupiter's  rotation, 

Io's  revolution,  and  the  earth's  rotation  respectively,  and  P„   is 

Sy 

the  Jovian  synodic  year.   For  j  =m=  1,  k=  -4,  and  n  =  0,  a  spurious 
periodicity  is  expected  at  13.43018  x  10   Hz ,  or  a  mean  synodic  period 
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of  20  40  59.2,  which  is  within  1.8  seconds  of  Bigg's  satellite  A. 
The  author  feels  that  this  could  also  account  for  the  prominent 
dip  in  probability  at  180°  in  Figure  49. 

One  must  conclude  from  this  study  that  there  is  very  little 
reason  for  believing  that  these  hypothetical  satellites  even  exist 

Amalthea 


A  search  for  a  possible  influence  on  the  radiation  by  Moon  V, 
or  Amalthea,  was  carried  out  by  Register  (1968)  and  the  author. 
The  method  of  analysis  was  similar  to  that  described  earlier  in  this 
chapter.   Data  at  the  frequencies  of  12,  18,  22.2,  27.6  MHz  were 
used  in  the  study  and  all  results  were  negative.   Bigg  (1969)  claims 
the  existence  of  a  small  effect  but  gives  no  data. 

From  these  studies  and  the  work  of  others  investigating  the 
remaining  three  Galilean  satellites  it  can  be  concluded  that  any 
effect  satellites  other  than  Io  may  have  on  the  emission  from  Jupiter 
must  be  very  small. 


CHAPTER  VII 
SUMMARY  AND  CONCLUSIONS 

The  pulse  character  analysis,  in  which  the  normal  and  millisecond 
pulses  were  independently  studied,  indicated  that  the  millisecond  pulses 
are  far  more  sensitive  to  the  position  of  Io  than  are  the  slower  pulses. 
It  also  showed  that  the  non-Io -related  component  of  source  A  is  composed 
almost  entirely  of  normal  pulses. 

The  study  of  the  millisecond  pulses  from  source  B  clearly 
showed  this  source  to  be  bifurcated,  and  it  also  indicated  that  quite 
different  positions  of  Io  are  required  for  the  stimulation  of  radiation 
from  the  two  subsources,  B..  and  B  . 

The  strong  influence  of  Io  on  the  millisecond  pulses  has  been 
used  by  the  Florida  group  in  preparing  prediction  tables.   These 
tables  have  been  very  useful  in  conducting  several  special  experiments 
in  which  only  the  millisecond  pulses  were  of  interest. 

No  firm  conclusion  can  be  drawn  concerning  the  frequency  dependence 
or  the  source  dependence  of  the  apparent  source  drift  because  of  the 
difficulty  in  locating  sources  B  and  C  at  the  lower  frequencies  and 
for  apparitions  of  low  activity. 

The  determination  of  the  rotation  period  of  System  III  using 
Duncan's  method  must  be  questioned  since  this  study,  using  essentially 
the  same  method,  derived  a  period  differing  from  Duncan's  by  0.33 
seconds,  a  discrepancy  far  greater  than  the  error  involved  in  the  method. 
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The  analysis  of  the  durations  of  the  storms  did  not  indicate 
any  systematic  variation  over  the  apparitions  and  thus  could  not 
account  for  the  observed  sinusoidal  drift  of  the  sources.   At  present, 
the  most  plausable  explanation  of  the  drift  is  the  one  due  to  Carr 
et_  al_.  (1970),  which  correlates  the  drift  with  the  declination  of  the 
earth  as  seen  from  Jupiter.   Since  the  period  of  the  drift  is  about 
12  years,  many  more  years  of  data  will  be  required  to  confirm  the 
theory  of  Carr  et_  al. 

In  the  final  question  considered,  that  of  possible  Io-type 
influence  on  the  radiation  by  other  satellites,  one  can  conclude  from 
this  study,  as  well  as  the  work  of  others,  that  any  such  influence 
must  be  very  small  in  comparison  to  that  of  Io . 


APPENDIX 
RECEPTION  PROBABILITY  HISTOGRAMS 
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